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We perform a comprehensive analysis of the spectrum of proper oscillations (quasinormal modes), 
transmission/reflection coefficients and Hawking radiation for a massive charged scalar field in the 
background of the Kerr-Newman black hole immersed in an asymptotically homogeneous magnetic 
field. There are two main effects: the Zeeman shift of the particle energy in the magnetic field 
and the difference of values of an electromagnetic potential between the horizon and infinity, i.e. 
the Faraday induction. We have shown that "turning on" the magnetic field induces a stronger 
energy-emission rate and leads to "recharging" of the black hole. Thus, a black hole immersed in 
a magnetic field evaporates much quicker, achieving thereby an extremal state in a shorter period 
of time. Quasinormal modes are moderately affected by the presence of a magnetic field which is 
assumed to be relatively small compared to the gravitational field of the black hole. 
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I. INTRODUCTION 

Magnetic field is one of the most important con- 
stituents of the cosmic space and one of the main sources 
of the dynamics of interacting matter in the Universe. 
Weak magnetic fields of about a few ijlG exist in galax- 
ies and clusters of galaxies, while very strong magnetic 
fields of up to 10 4 — 10 8 G are supposed to exit near su- 
permassive black holes in the active galactic nuclei and 
even around stellar mass black holes [H-ll] . Magnetic field 
near a black hole leads to a number of processes, such as 
extraction of rotational energy from a black hole, known 
as the Blandford-Znajek effect Q , the charging of a black 
hole due to accretion of charged matter yjj, the forma- 
tion of an induced electric field on the black hole surface 
, negative absorption (masers) of electrons Q , and so 
on. In addition to stellar mass and galactic black holes, 
miniature black holes could be immersed in a strong mag- 
netic field if created in a laboratory or observed in cosmic 
showers 

Even a relatively weak magnetic field can considerably 
affect the behavior of charged particles/fields due to usu- 
ally not weak coupling eB between the particle charge 
e and the magnetic field B. Therefore, charged massive 
fields are interesting models for theoretical study of inter- 
action of a magnetized black hole with its surroundings. 
As the simplest case, one may neglect the spin of the field 
and consider the complex massive charged scalar field. 
Still, the interaction of particles due to the spin can also 
significantly affect the particles' state and thus deserves 
a separate consideration. As a first step in this direction, 
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we shall consider here a charged rotating black hole, given 
by the Kerr-Newman solution, and a charged massive 
scalar field propagating on its background and immersed 
in an asymptotically homogeneous magnetic field. The 
magnetic field is supposed to be weak enough so that the 
metric does not deviate from the Kerr-Newman one, i.e. 
the magnetic field does not distort the geometry of the 
space-time but only interacts with other electromagnetic 
charges in the system. 

Particles and fields in the vicinity of a black hole 
slightly change the background space-time of a system. 
Therefore, the addition of a field to a black hole space- 
time can be considered as a perturbation. At the clas- 
sical level, the perturbation can be described by its 
damped characteristic modes, called the quasinormal 
modes (QNMs) [H], EH and by the scattering properties, 
which are encoded in the S-matrix of the perturbation. 
Quasinormal modes are proper oscillations of the pertur- 
bation which dominate at late time in the response of a 
black hole to z perturbation. The complex frequencies of 
such oscillations do not depend on the manner of exci- 
tation but only on the parameters of the black hole and 
the field under consideration. Therefore, they are usually 
called the "fingerprints" of a black hole. 

In the same way, as quasinormal modes are an essen- 
tial classical characteristic of a black hole, the thermal 
Hawking radiation is its essential quantum feature that 
carries information about the dynamics of evaporation of 
the black hole. For large astrophysical black holes, the ef- 
fect of Hawking evaporation is certainly negligible for the 
black hole dynamics but not for the behavior of particles 
in its vicinity. Emission of Hawking radiation is signif- 
icant for primordial black holes and huge for miniature 
black holes which are considered in the higher dimen- 
sional gravity and string theory. According to brane- 
world scenarios, our world is assumed to be a (3+1)- 
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dimensional brane which is embedded in a higher dimen- 
sional bulk. Gravity is supposed to be much stronger 
at small distance creating possibilities for the formation 
of an event horizon in particle collisions even at energies 
achievable at the Large Hadron Collider |-sl . Estimates 
show that, once such a black hole is created, it will almost 
immediately evaporate, so that the life-time of the minia- 
ture black holes is about 10~ 25 — 10~ 32 sec. Although 
created in this way miniature black holes would be intrin- 
sically higher dimensional. Nevertheless, our considera- 
tion here of the 4-dimensional black holes can be useful 
for two reasons: First, it is known that the higher dimen- 
sional black hole emits particles mainly "on the brane," 
i.e. the process of quantum evaporation of the higher 
dimensional black holes is probably qualitatively simi- 
lar in many aspects to the one for 4-dimensional black 
holes. Second and more important, in the proposed ap- 
proach, we can realize how a magnetic field can influence 
the process of evaporation. Thus, when talking about 
QNMs and Hawking radiation of black holes immersed 
in a magnetic field we have in mind not only astrophys- 
ical black holes but also primordial and miniature black 
holes. 

Thus, we could say that the quasinormal spectrum and 
Hawking radiation are, respectively, classical and quan- 
tum "fingerprints" of a black hole. Quasinormal modes 
and Hawking radiation also have one technical point in 
common: analysis of QNMs as well as of the Hawking 
radiation (in semiclassical approximation) begins from 
the linear perturbations of the fields under consideration 
whose dynamics should be reduced to a single wavelike 
equation, called the master equation. 

Up to now there are two kinds of analysis of quasinor- 
mal modes which, in a sense, are complimentary to the 
present work. Quasinormal modes of the massive and 
massless charged scalar field around charged black holes 
(without a magnetic field) were studied in (16l - |l8j . Quasi- 
normal modes of a neutral scalar field around black holes 
immersed in a strong magnetic field were calculated in 
[H [H. In [H HI, the black hole was described by the 
Ernst-Schwarzschild solution which contains a magnetic 
field as a parameter because the magnetic field is implied 
to be strong enough in order to deform the black hole 
geometry significantly. However, such strong geometry- 
deforming magnetic fields have little probability of exist- 
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ing in nature Q. 

Here, we shall consider a more realistic situation and 
assume that the magnetic field is not strong enough to 
deform the Kerr-Newman black hole metric. The correla- 
tion of the quasinormal frequencies, the reflection coeffi- 
cients and the energy emission rates with the parameters 
of the black hole (mass M, charge Q, angular momen- 
tum a) and of the scalar field (mass fj,, charge e) are 
analyzed here through a comprehensive numerical study. 
The Hawking radiation for charged particles (without a 
magnetic field) around nonrotating and rotating black 
holes was considered in [To| for 4-dimensional black 
holes and in [n], [l2| for higher dimensional scenarios. 

In the system under consideration, the coupling of par- 
ticle charge e with the magnetic field B leads to the 
Zeeman shift of the energy /i 2 —> /i 2 — eBm [HI]. The 
rotation of the black hole in the magnetic field, in its 
turn, leads to the appearance of the induced charge on 
the black hole surface and to the difference in values of 
the electromagnetic potential at the horizon and at in- 
finity, that is, the Faraday induction. We shall observe 
how these two effects, the Zeeman effect and the Faraday 
induction, are reflected in the processes of classical and 
quantum radiation. Qualitatively, these two effects were 
considered in the vicinity of the Kerr-Newman black hole 
by Galtsov and collaborators [6j. Here we shall give an 
accurate quantitative analysis for the above case. We 
shall calculate characteristic quasinormal modes, reflec- 
tion/transmission coefficients, and the emission rates for 
Hawking radiation of the charged massive scalar field in 
the background of the Kerr-Newman black holes and in 
the vicinity of the asymptotically homogeneous magnetic 
field. 

The paper is organized as follows. Section II is devoted 
to the separation of variables for the scalar field in the 
Kerr-Newman background under nonzero magnetic field. 
In Section III we describe the numerical procedures for 
finding eigenvalues of the separated angular equation and 
quasinormal frequencies and reflection/transmission co- 
efficients. Section IV is devoted to calculations of the 
quasinormal modes. Section V discusses classical scat- 
tering and calculates the energy and momentum emis- 
sion rates for the Hawking radiation. We present our 
conclusions in Sec. VI. 



II. WAVE-LIKE EQUATION 

The Kerr-Newman metric can be written in the following form: 

ds 2 = ^(dt -asm 2 dd(f>) 2 - ^-^-[adt - {r 2 + a 2 )d(f>} 2 - ^-dr 2 - Ed(9 2 , (1) 
Zj Zj ZA 

where 



A = r 2 -2Mr + a 2 + Q 2 , E = r 2 + a 2 cos 2 9. (2) 
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Here M is the black hole mass, Q is its charge, and a is the angular momentum per unit mass. The event horizons 
are situated at 



r = r ± =M± y/M 2 -a 2 -Q 2 . (3) 

In the above description we have not taken into account the influence of the magnetic field onto the black hole 
background. Under these conditions the background electromagnetic field can be written as 

A = A^dx* = —{dt- a sin 2 Od<j>). (4) 

The KN metric does not depend on the coordinates t and <f>, so that there exist the two Killing vectors £m = (1, 0, 0, 0) 
and £(0) = (0, 0, 0, 1). One can see that the Killing vectors for vacuum metrics satisfy the same equations as the 4- 
potentials A^. This suggests the following form of the 4-potential: 



A" = -B 
2 



ft) + 2 <) 



2M$r (5) 



The gauge transformations, 



d ( (Q-1aMB\ \ 

reduce the 4-potential to the Coulomb form A$(r = oo) = 0. Thus, the "full" electromagnetic 4-potential of the 
system which includes the electric field of the charged black hole as source and a magnetic field B living "in the 
background" of a Kerr-Newman black hole has the form 

_ f (Q-2aMB)r(r 2 +a 2 ) -B (Q - 2aMB)ra 
A ~\ AS » »°>2-+ £x 

„ „„ B 2 sm 2 9r /9 9n9 a , 2nl ^B , „ (Q - 2aMB) 2 r 2 
A ll A»= — — [{r 2 +a 2 ) 2 - Aa 2 sm 2 9] - (Q - 2aMB) —r sin 6 + — — -!■ . (7) 

The Klein-Gordon equation for the charged massive scalar field in the vicinity of a Kerr-Newman black hole and in 
the presence of the homogenous magnetic field B has the following general covariant form [l3j : 

0""(V M + *eA M )(V„ + ieA*)* + ^ = (8) 
Since V^A' 1 = 0, the latter can be reduced to the following form [l3T |: 

8 / A 3#\ 1 8 ( f(r 2 +a 2 ) 2 , , \ 9 2 * / r 2 + a 2 \ 8 2 ^ ( 1 a 2 \ 8 2 ^ 
A— H -— sin6>— - : — a 1 sin 2 9 — r + 2a 1 - 1 ' 1 



9r V dr / sind 86 \ 89 J \ A ) 8t 2 \ A / di<90 Vsin 2 <9 A/ <9 2 </> 

>(r 2 +a 2 )(Q - 2aMB) 8^ ({Q-2aMB)ra BT,\ 



2ie 



A 9t V A 2 



+ (e A^A^ - ^ 2 )E* = 0. (9) 



As was shown in [ijj, the separation of radial and angular variables in the whole space is impossible for this 
equation. Nevertheless, if one considers only the region which begins at the event horizon and ends at some distance 
far from the black hole r » r +1 and uses the following approximations: 

eBr 2 + < 1, eQ < 1, (10) 

then we get, 

= (Q - 2aMff) 2 r 2 
M AE 

Under these conditions the separation of variables is possible, that is, we can write 



»*+™<l>S(6)R(r)/y/r 2 + a?. (11) 
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Here the function S(9) obeys the equation 

/ d 2 d m 2 \ 

\^+ cot9 ™ 2 ( au} f sin 2 9 + 2mauj + A - (> 2 - eBm)a 2 cos 2 9 S{9) = 0. (12) 

\oa z oO sin 9 J 

This equation can be solved numerically for any value of lu in the same way as the equation for massive scalar field 
in the Kerr black hole background [2l| with the effective mass /i 2 ff = /i 2 — eBm. We note that, when fi e g — 0, Eq. 
(Ti"2"]) reduced to the well-known equation for the spheroidal functions. In this case, the separation constant \(u>) can 
be found numerically using the continued fraction method [22]. When the effective mass is not zero, the separation 
constant can be expressed, in terms of the separation constant for spheroidal functions, as 

\(lu, /i G ff ) = A(f2) + 2ma(tt — lu) + /i 2 ffa 2 , 

where 



fl = W lu 2 — fi 2 s = \J 'lu 2 — /i 2 + eBm. 

The sign of Q might be chosen here arbitrarily, but we fix it so that Re(fl) and Re(to) are of the same sign. This 
allows one to recover easily the limit of /i e g = and, later, simplify fixing of the boundary conditions for the radial 

part. When a — 0, one can find that A = 1(1 +1), I = 0, 1, 2 For the nonzero values of a, the separation constant 

can be enumerated by the integer multipole number I > |m|. 

Using the new tortoise coordinate r* , the radial part can be written as a wavclike equation, 

d?_ 

^dr 2 

where 



V(r) R(r) = 0, (13) 



A 



and where the effective potential has the following form 



A / „ , , , (rAY 3Ar 2 \ ( ma Q - 2aM B\ 2 , . 

V ( r ) = ( 2, 2N2 A + M 2 - eBm)r 2 + - - « - - T — - - er ^ . (14) 

From this form of the effective potential one can realize the two effects: the mass of the field gained an effective 
term fi 2 — eBm, and the black hole charge gained an addition as well, Q — >• Q — 2aMB. The first effect is the 
well-known Zeeman effect, which is the shift of energy of a charged particle (with a charge e) in the magnetic field 
due to interaction of a magnetic field B with an azimuthal momentum m. In systems which are more symmetric than 
ours, i.e. with degenerated m-states, the Zeeman effect leads to splitting of the m-degeneration, and is well known in 
quantum mechanics. In the case when the effective potential allows for nondegenerated m-states, the Zeeman effect 
simply corresponds to a shift in the particle's energy. 

The second effect is more remarkable. Once a rotating black hole is immersed into a magnetic field, the electrostatic 
potential between the horizon and infinity acquires a difference due to the presence of a magnetic field, which is 

5 A — Ah A r - — 2aMB 

In other words the black hole receives an additional induced electrostatic force Find = 2aMB jr 2 + . This is nothing but 
the Faraday induction. It should be noted that this effect can be significant even for neutral black holes Q = 0, and 
can be applied to large astrophysical black holes which cannot possess large electric charge. 
The asymptotics of the effective potential near the horizon and at infinity are 



V(r) -> -n 2 , r — v oo, SI = y/co 2 - p 2 + eBm, (15) 

_ 2 . ma + er + (Q - 2aMB) 

Vvr) — > — to , r — >■ r+, to = to 5 5 . (lb) 

a z + r+ 

During super-radiance the black hole can be charged until it reaches the "extremal" charge Q = 2aMB. 

In the following sections we will study the influence of the above-mentioned Zeeman and Faraday effects on the 
classical (quasinormal) and quantum (Hawking) radiation of black holes. However, before starting the numerical 



5 



study of the wavelike equation (|13p . let us mention one more constraint related to our analysis. If the electric field 
2aMB/r + induced on the horizon is as strong as the Schwinger field /x 2 /e, then the electrodynamic process of particle 
production is initiated and continues until the maximum value of the charge is reached. This maximum values of the 
charge is Q — 2aMB. Indeed, imagine vacuum as consisting of virtual pairs e + e _ where electrons and positrons, after 
transforming to the real ones, become separated by distance of the order of the Compton wavelength A = 2ir[i~ 1 . 
If the work done by the electric field, which is eEX, is as large as the rest mass of the two particles 2/x, then the 
virtual pair turns into a real one: 2n[i~ 1 eE > 2/i and consequently E ss /i 2 /e, where E is the induced electric field. 
We did not take this Schwinger mechanism into consideration in the present paper. However, in the conclusion, we 
shall suggest simple arguments showing that the Schwinger mechanism of pair production will enhance the process of 
"recharging" a black hole and force the black hole to evaporate faster. 



III. NUMERICAL METHODS 

In this section we shall briefly discuss the two clas- 
sical numerical methods (Frobenius and WKB) used 
for calculations of the quasinormal modes and the 
shooting method used for calculations of the transmis- 
sion/reflection coefficients. 



The function y{r) must be regular at the horizon and 
spatial infinity and 



\Jr 2 + a 2 x - / r - /■_ 

y ( r ) = }^ a k 



fe=0 



The coefficients satisfy the three-term recurrence re- 
lation. 



A. Quasinormal modes 



a„a 



n^n+l 



+ P„a n + ~f n a n - 1 = 0, n > 0, 7o = 0, (17) 



In order to calculate quasinormal modes we impose the 
quasinormal mode boundary conditions for the wave Eq. 
([T3| . i.e. we require that at the black hole horizon we 
have only purely ingoing waves, 



-oo) cx exp(— iuir*), 



while we should have only purely outgoing waves at spa- 
tial infinity, i.e. 

i?(r* — > oo) cx exp(if2r*). 

Thus, no waves are coming from the horizon or infinity, 
which implies that w are proper oscillation modes in the 
black hole response to an "instantaneous" perturbation. 
In other words, when the perturbation decays, the source 
of the initial perturbation is not acting anymore. 

Equation (|13[) has an irregular singularity at spatial 
infinity and four regular singularities at r = r+, r = 
r_ = (Q 2 + a 2 )/r + and r — ±ia. The four singularities 



appear due to the prefactor (r 2 



-1/2 



in CLE]). The 



appropriate Frobenius series is determined as 



R(r) 



r+ 



-iui I AttTh 



e iQr (r-r_)S(r), 



where 



(7= 



/i 2 — eB(m — 2aw) \ 
2Q ) 

and Th is the Hawking temperature 

A'(r + ) 



(r+ +r_), 



T H = 



47r(r 2 _ + a 2 ) ' 



where a n , [3 n , 7„ can be found in an analytic form. We 
do not write these coefficients here because they have 
quite a cumbersome form. Notice that the factor ^-t^ 
removes the singularities r — ±ia of y(r). Except for 
this factor, we would have had a five-terms recurrence 
relation due to the additional singular points. 
By comparing the ratio of the series coefficients 



a n 
dn+l 



7« 



Oin-l 



&n Pn- 



/?„_2-a„_37„_ 2 /.- 
7n+l 



0n 

On 



n+1 



Qn+l7n+2 
/3n+2— "n+27n + 3/- 



(18) 



we obtain an equation with a convergent infinite contin- 
ued fraction on its right side: 



(3n 



-lln 



(19) 



0n—2 — Q!»-37r(-2/... 



n+1 



Qn + l7Ti + 2 



which can be solved numerically by minimizing the ab- 
solute value of the difference between its left and right 
sides. Equation ([T9f has an infinite number of roots, but 
the most stable root depends on n. Generally the larger 
number n corresponds to the larger imaginary part of the 
root w [HI . 

Note that the case under consideration allows one to 
use the Nollert procedure (24[, in order to improve con- 
vergence of the infinite continued fraction, which is useful 
for searching roots with a very large imaginary part. 

For an additional check of the accurate numerical re- 
sults obtained by the convergent Frobenius method, we 
shall use also the WKB formula of the 6th order beyond 
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the eikonal approximation [25l . \2t 
following form: 



The formula has the 




n = 0,1,2. 



(20) 



and the correction terms A$ were obtained in [25|, |26j and 
depend on higher derivatives of V at its maximum with 
respect to the tortoise coordinate r*, and n labels the 
overtones. The WKB approach was developed by Schutz 
and Will HH and extended to the 3rd [25j and 6th HI 
orders. It can be effectively used not only for finding low- 
lying quasinormal modes (see, for instance, [13, [29| and 
references therein), but also for calculations of the trans- 
mission/reflection coefficients in various problems j30|- 



B. Reflection coefficients 

For calculations of the emission rates of particles due to 
Hawking radiation, one needs first to solve the problem of 
classical scattering in order to obtain the gray-body fac- 
tors. This implies the posing of classical scattering bound- 
ary conditions. At the event horizon, this again means 
imposing the boundary condition which corresponds to 
a purely ingoing wave, while at spatial infinity (r — > oo) 
we have a different condition from the one used for the 
quasinormal modes, 

R(r) ~ Z in cxp(-if2r*) + Z out exp(«f2r*), 

where Zi n and Z out are integration constants which cor- 
respond to the ingoing and outgoing waves respectively. 
Thus, we would like to know which portion of particles 
will be able to pass through the barrier of the effective 
potential. 

Introducing the new function 



P(r) = R(r) 



iuj /AttTh 



and choosing the integration constant as P(r + ) = 1, we 
expand Eq. (|13l) near the event horizon and find P'(r + ), 
which completely fixes the initial conditions for the nu- 
merical integration. Then, we integrate Eq. (1131) numer- 
ically from the event horizon r+ to some distant point 
Tf 3> r + and find a fit for the numerical solution far 
from the black hole in the following form: 



P(r) = Z m P in (r) + Z out P out (r), 



(21) 



where the asymptotic expansions for the corresponding 
functions are found by expanding (Tl3|) at large r as 



Pn{r) 



e -mr r -ia ( 1 + p 0-) r -i +p (?) r -2 



Pout (r) = ^^(l + Pii^+Pflr- 2 + ...). 

The fitting procedure allows us to find the coefficients Zi n 
and Z out . In order to check the accuracy of the calculated 




FIG. 1: The effective potentials for the charged scalar field 
(I = m — 2) in the Schwarzschild background (r+ = 1) with 
the magnetic field eB = 0.04 (top blue curve), eB — 0.1 
(middle green curve), eB = 0.25 (bottom red curve). One can 
see that the potential is negative only deeply in the region 
r 3> 1/V eB, which is beyond the region of validity of the 
approximation (fl0|l . 



coefficients, one should increase the internal precision of 
the numerical integration procedure, the value of rf, and 
the number of terms in the series expansion for Pi n (r) 
and Pout (f) ) making sure that the values of Zi n and Z out 
do not change within desired precision. 

If the coefficients Zi n and Z out are calculated, one can 
find the absorbtion probability 



1 — \Z ou t/Zi, 



(22) 



This will be used later for calculations of the emission 
rates for energy momentum and charge of the black hole. 
This approach was also used for analysis of Hawking ra- 
diation of higher dimensional simply rotating black holes 
[32| , and of Gauss-Bonnet black holes [33J , and it showed 
an excellent agreement with the analytical approach. 



IV. QUASINORMAL MODES 

First, let us briefly review previous works on quasinor- 
mal modes of Kerr-Newman black holes. Apparently the 
first work on QNMs of Kerr-Newman black holes was by 
one of us [3(| , where gravitational perturbations with a 
frozen Maxwell field were considered. More accurate nu- 
merical results for gravitational perturbations were pre- 
sented in 37]. Quasinormal modes of the Dirac field 
were considered in 38]. The case of charged scalar and 



Dirac modes was analyzed in the background of Kerr- 
Newman black holes allowing for a positive cosmological 
constant [lH . All the above papers observed no unstable 
modes in the quasinormal spectrum (except the superra- 
diant ones) , although unstable modes of a charged scalar 
field were found for an asymptotically anti-dc Sitter back- 
ground Hi]. 

Before we examine the dependence of QNMs on vari- 
ous parameters of the system, we should first discuss the 
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TABLE I: QNMs of the massless scalar field in the background 
of a nonrotating uncharged black hole, eB = 0.05. The Frobc- 
nius method gives the unstable mode, which does not appear 
when we use the WKB method, supposing an asymptotically 
flat background. 



mode 


unstable 


stable 


WKB fundamental 


£ = m = l 
£ = m = 2 
£ = m = 3 
£ = m = 4 
£ = m = 5 


0.2236i 
0.3162i 
0.3873i 
0.4472i 
0.5000i 


0.5747 - 0.2020i 
0.9516 - 0.1988i 
1.3331 - 0.1973i 
1.7162 - 0.1963i 
2.0999 - 0.1956i 


0.5747 - 0.2022i 
0.9515 - 0.1989i 
1.3331 -0.1973? 
1.7162 - 0.1963i 
2.0999 - 0.1956i 



Im(o>) 




FIG. 2: Real (left panel, graphs from top to bottom) and 
imaginary (right panel, graphs from bottom to top) parts of 
the fundamental (n = 0) QNM as a function of eB for £ — 
m = 0, Q = 0, a=0.2 (blue), a=0.4 (green), a=0.6 (orange), 
a=0.8 (red), a=0.99 (magenta). 



stability of the system. The effective potential Eq. (fT4|) 
contains a term proportional to /i 2 — eBm which works 
as an effective mass term, and, when /i 2 < eBm, this 
term is negative. It is well known that a massive scalar 
field with negative /j 2 is unstable even for tiny negative 
values of the square of mass (34|. Thus, the instabil- 
ity is expected when formally considering exact solutions 
of the wave equations (|13j) . but certainly not for a real 
physical situation. The reason is that the instability due 
to negative square of mass comes from infinite negative 
fall-off of the effective potential for scalar field at spa- 
tial infinity. In our case, however, "infinity" is located at 
r ~ r + (eB)^ 1 , and further from this distance the wave 
equation is not valid because of the approximation (|10p 
which has been used for the separation of variables. In 
Fig. [U one can see that the effective potential is positive 
definite in the region of its validity and is negative only 
for values of r which are seemingly larger than (eB) -1 . 

There are reasons to expect that the true effective po- 
tential for the black hole immersed in an asymptotically 
uniform magnetic field will inevitably lead to instability 
due to the infinite energy of the magnetic field. Anal- 
ysis of particle motion around Ernst-Schwarzschild and 
Ernst-Kerr black holes shows that the effective potential 
for such particles diverges at infinity. This means that 
the magnetic field which fills in all the Universe will cre- 
ate an effective confining box. Thus, at infinity it will 
be appropriate to use Dirichlet boundary conditions. A 
rotating black hole in such a confining box will inevitably 
be unstable through the mechanism of superradiance. In 
a real world the magnetic field is certainly assumed to 
vanish at infinity, so that no confining box will appear. 
When using approximation (fTU|) . we "cut" the effect of 
the confining box at infinity in a natural way. 

In nature, infinity means a region far from the black 
hole r r + which can be approximately treated as 
asymptotically flat. In practice, one should match the 
considered solution with asymptotically homogeneous 
magnetic field at "infinity" r rs r + {eB)~ l with some 
asymptotically flat solution. Fortunately, as was shown 
in [35j, the major scattering properties of fields, including 
the low-lying quasinormal modes, depend on the behav- 
ior of the effective potential only in some region near 
the black hole (if the black hole is not an anti-de Sitter) , 
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FIG. 3: Real (left panel, graphs from top to bottom) and 
imaginary (right panel, graphs from bottom to top) parts of 
the fundamental (n = 0) QNM as a function of eB for £ = 1, 
m — 0, Q — 0, a=0.2 (blue), a=0.4 (green), a=0.6 (orange), 
a=0.8 (red), a=0.99 (magenta). 




FIG. 4: Real (left panel, graphs from top to bottom, except 
for a — 0.2) and imaginary (right panel, graphs from bottom 
to top) parts of the fundamental (n = 0) QNM as a function 
of eB for I = 1, m = 1, Q = 0, a=0.2 (blue), a=0.4 (green), 
a=0.6 (orange), a=0.8 (red), a=0.99 (magenta). 
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FIG. 5: Real (left panel, graphs from top to bottom) and 
imaginary (right panel, graphs from bottom to top) parts of 
the fundamental (n = 0) QNM as a function of eB for I = 1, 
m = —1, Q — 0, a=0.2 (blue), a=0.4 (green), a=0.6 (orange), 
a=0.8 (red), a=0.99 (magenta). 
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while the form of the effective potential far from the black 
hole has no impact on the results. This can easily be ex- 
plained because the process of scattering occurs mainly 
near the maximum of the potential barrier. 

When computing QNMs with the help of the Frobe- 
nius method, we do not take into account this peculiar- 
ity of infinity and treat infinity as a mathematical one. 
Thus, in addition to a number of damped stable modes, 
we must find some "unstable" modes by the Frobenius 
method. Indeed, we find and tabulate them in Table I, 
where one can also find estimations of QNMs derived us- 
ing the WKB method. Unlike the Frobenius method, the 
WKB formula 26] implies that one has a positive defi- 
nite decaying potential at infinity. Therefore, the WKB 
formula mainly approximates the behavior near the max- 
imum of the potential and does not reproduce those "un- 
stable" modes. That is physically adequate as the grow- 
ing modes appear within Frobenius approach only due to 
improper "extension" of the wave equation (| 13[) outside 
the region of its validity. 

Finally, another instability occurs due to the so-called 
superradiance: the massive field has a local minimum far 
from the black hole which works as an effective potential 
wall, so that the wave amplified due to extraction of ro- 
tational energy of the black hole (superradiance) can be 
reflected back from the distant wall. Repetition of this 
process leads to unbounded growth of the perturbation. 
Superradiant instability is shown to be always negligi- 
bly small [29j , so that the evaporation time of miniature 
black holes is much shorter than the characteristic time 
of the instability growth. For large, astrophysical black 
holes, superradiant instability of massive fields means 
that the quantum field will go over to the higher non- 
superradiant state. In addition, unstable modes are not 
in the quasinormal sector of the black hole spectrum. 
Therefore, we do not need to give a detailed analysis 
of superradiant modes here. Moreover, such an analysis 
would be technically inaccurate within our approach be- 
cause the effective potential is known only in some prox- 
imity of the black hole and is not exactly known far from 
the black hole where the local minimum is localized. 

Quasinormal modes of the Kerr-Newman black hole 
immersed in a magnetic field for massive charged scalar 
field will be determined by a number of parameters, seven 
altogether: the black hole parameters Q, M, and a, the 
magnetic field B, the scalar field parameters /1 and e, 
and its quantum numbers m and £. Therefore, com- 
plete investigation of the quasinormal modes correlation 
on these parameters would include an enormous amount 
of numerical data. We shall show here only the most 
representative plots for dependence of QNMs on various 
parameters. We present all our quantities in units of the 
black hole horizon. 

In Figs. [3J [3J HI [5] one can see for the Q = case 
that the m = modes and the modes with nonvanish- 
ing azimuthal number behave quite differently. Actually, 
modes with m = (Figs. [21 EJ) have decreasing damping 
rate as the angular momentum per unit mass a increases. 



In the regime of relatively small values of eB the damping 
rate increases roughly linearly with eB. The real oscilla- 
tion frequency Re{uj) decreases linearly with eB and also 
decreases for growing a. We can also see that eB coupling 
has greater influence on Re(ui) than on Im(ui), which re- 
mains almost unchanged within the region of small eB. 

For modes with m > 0, both Re(cu) and Im(ui) lin- 
early decrease with eB (Fig. [4]), with one peculiarity: for 
moderate negative values of eB the Re(uj) is not mono- 
tonically decreasing with a for all eB anymore. In a 
large region of both positive and negative values of eB, 
the Re(ui) monotonically decreases with a up to some 
minimal value and then increases. This explains the in- 
tersection of curves in Fig. 21 For negative m, we did 
not observe such a minimum (Fig. [5]) and the behavior 
is quite similar to the m — case. 

In Fig. [HI one can see that Re(uj) of the fundamental 
mode (£ = in — 0) is monotonically decreasing with Q 
for eB < and monotonically increasing for eB > 0. For 
eB < 0, Im(uj) monotonically grows with Q, while for 
eB > 0, Im{uj) as a function of charge Q decreases up to 
a minimum at some moderate value of Q, and then starts 
growing. As in the limit m = and Q = 2aMB Eqs. (fT2")) 
and (I13[) do not depend on the field charge e, the upper 
and lower curves in Fig. [5] coincide for Q — 2aMB. In 
Fig. [6l the quasinormal behavior for £ = l,m = is 
similar to Fig. [71 where modes with positive and neg- 
ative e also coincide in the limit Q = 2aMB. Modes 
with I = m = 1 have also monotonically decreasing (in- 
creasing) Re(uj) as a function of the charge Q for eB < 
(eB > 0), while Im(uj) is monotonically increasing for 
both positive and negative eB (see Fig. [7J. The same 
monotonic growth of Im(uj) happens for I = l,m = — 1 
mode, so the Re(uj) has an opposite behavior: it grows 
for eB < and decreases for eB > 0. 

Finally, let us discuss the dependence of quasinormal 
modes on the charge of the field e. In Fig. [8j one can 
see £ = 0, 1 m = modes and £ = 1, m = ±1 modes 
with a charge Q which is equal to the "extremal" value 
2aMB. Modes with m — naturally form an almost 
horizontal line because e enters into the wave equation 
in combination with m or Q — 2aMB. Thus, exactly in 
the limit Q = 2aMB, we have a single mode which is 
independent of e. For I = 1, m = 1, i?e(w) decreases 
as a function of e, so this decrease has some small local 
peaks at larger negative values e. In a similar way, 1=1, 
m = — 1 modes have both a real and an imaginary part of 
lo which almost monotonically increase up to small peaks 
for moderate values of |e|. For sufficiently small values 
of |e|, the dependence of w on e is strictly monotonic. 

Summarizing, we can say that the Zeeman effect and 
the Faraday induction influence the quasinormal spec- 
trum in a rather complicated way, where one cannot eas- 
ily distinguish these two effects. This happens because 
the magnetic field acts in a different way on modes with 
different azimuthal numbers m: m = modes are usu- 
ally only slightly affected by the magnetic field, while 
m 7^ are strongly influenced due to the extra coupling 
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FIG. 6: Real and imaginary parts of the fundamental (n = 0) QNM as a function of Q for Q = 0, 2aMB/3, aMB, 2aMB, 
B = 0.2, a — 0.6, jj, = 0, e = +0.2: £ = (red, right panel: bottom), I = 1 (magenta, left panel: top) and e = —0.2: £ = 
(blue, left panel: bottom), 1=1 (green, right panel: top), m = 0. 
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FIG. 7: Real (left panel) and imaginary (right panel) parts of the fundamental (n = 0) QNM as a function of Q for Q = 
0, 2aMB/3, aMB, 2aMB, B = 0.2, a = 0.6, fj, = 0, from top to bottom: < = m = le = +0.2 (red) and e = -0.2 (blue), 
t = — m = 1 e = —0.2 (green) and e = +0.2 (magenta). 
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FIG. 8: Real (left panel) and imaginary (right panel) parts of 
the fundamental (n = 0) QNM as a function of e for fi = 0.1, 
a = 0.25 Q = 0.9, B = 1440/749 w 1.92, J = m = (blue, 
bottom horizontal) 1 = 1, m = (red, top horizontal), Z = 
1, m = 1 (cyan, top), 1 = 1, m = — 1 (green). 



term Bm. Also, the quasinormal behavior is different 
for negative and positive charge of the field e. One com- 
mon feature of quasinormal modes is that the larger the 
magnetic field B is, the longer the QNMs live, if m < 0, 



while for m > 0, they live less long. A similar longer life 
of the quasinormal modes (yet, for all m and not only for 
m < 0) has been observed for the Ernst solution, i.e. for 
black holes immersed in a strong ma gne tic field that de- 
forms the geometry However, in [l9j a neutral scalar 
field was considered, so that the mechanism that influ- 
enced the quasinormal spectrum was the induced defor- 
mation of the space-time geometry due to magnetic field 
and not to the Zeeman shift and the Faraday induction. 
For the Ernst black hole a superradiant instability can 
also be considerably enhanced [2(| for huge values of the 
magnetic field strength, which seems to be very difficult 
to achieve. Quasinormal modes of charged scalar fields 
for nonrotating black holes have also been considered in 
a number of papers [l6l - IT8| . 
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FIG. 9: Grey-body factor of the Kerr-Newman black hole (left panel: a = 0.5, Q = 0.1, right panel: a = 0.5, Q = 0.5) in 
the magnetic field B = 2/3 due to massless charged particles (e = 3/20, £ = 1). From left to right: m = —1 negative charge 
particles, m = — 1 positive charge particles, m = positive charge particles, m = negative charge particles, m = 1 positive 
charge particles, m = 1 negative charge particles. For the "extremal" black hole m = grey-body factors for the particles and 
antiparticles are the same (black line). 
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FIG. 10: Energy-emission rate of the Kerr-Newman black hole (left panels: a = 0.5, Q = 0.1; right panels: a — 0.5, Q = 0.5) 
without the magnetic field (top panels) and with the magnetic field B — 2/3 (bottom panels) due to massless charged particles 
(e = 3/20). The red (top) and blue (bottom) lines correspond respectively to the same and the opposite signs of the charges 
of the black hole and the emitted particles. 
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FIG. 11: Angular momentum emission rate of the Kerr-Newman black hole (left panels: a = 0.5, Q = 0.1; right panels: a = 0.5, 
Q = 0.5) without the magnetic field (top panels) and with the magnetic field B = 2/3 (bottom panels) due to massless charged 
particles (e = 3/20). The red (top) and blue (bottom) lines correspond respectively to the same and the opposite signs of the 
charges of the black hole and the emitted particles. 



V. SCATTERING AND HAWKING RADIATION The energy-emission rate is 



A classical black hole in equilibrium does not emit any- 
thing. Nevertheless, when considering quantized fields 
around the black hole, the Hawking radiation appears: a 
black hole can create pairs of particles from the vacuum 
on the edge of its horizon. Particles with negative en- 
ergy go beyond the horizon, while particles with positive 
energy partially leave the black hole. When analyzing 
Hawking radiation of black holes we shall assume that the 
black hole is in thermal equilibrium with its surroundings 
in the following sense: the black hole temperature does 
not change between the emission of two consequent par- 
ticles. This implies the canonical ensemble as a model 
for the system. 

Not all positive energy particles can leave the back 
hole: part of them is reflected from the potential barrier 
surrounding the black hole. Thus, the energy-emission 
rate depends on the grey-body factors which give the 
fraction of particles penetrating the barrier. 

The emission rates for the energy, charge and angu- 
lar momentum are proportional to the grey-body factors. 



dM ~ ' y 

--*=£ e yiA 

the charge-emission rate is 



exp(w/T ff ) - 1 2tT 



(23) 



dQ 
dt 



oo £ 

EE 

g = m=-l 



\At, r 



duj 



(24) 



exp(Q/T H ) - 1 2vr 
and the angular-momentum emission rate has the form 
dJ V--V f , . ,9 fn du> 



dt 



E E / ia 



exp(w/T ff ) - 1 2tt" 



(25) 



Here, we perform the summation over all the possible 
values of the quantum numbers I and m. The grey-body 
factors are shown in Fig. |9]as functions of u>. There one 
can see that, for negative to, grey-body factors are larger 
for negatively charged particles than for positively ones, 
while, for m > 0, on the contrary, positively charged 
particles have larger grey-body factors than negatively 
charged ones. This is indirect influence of the Zeeman 
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FIG. 12: Charge-emission rate of the Kerr-Newman black hole (left panels: a = 0.5, Q = 0.1; right panels: a = 0.5, Q = 0.5) 
without the magnetic field (top panels) and with the magnetic field B — 2/3 (bottom panels) due to massless charged particles 
(e = 3/20). 



term eBm whose contribution depends on the sign of m 
and e. At first glance, negatively (relatively, the black 
hole charge) charged particles which are emitted radially 
should have smaller transmission coefficient than posi- 
tively charged ones: electromagnetic attraction of oppo- 
site charges diminishes the transmission of negative par- 
ticles. For particles which are radiated in all possible 
direction, this is certainly not so strict and the coupling 
with the azimuthal number m becomes important. 

Let us first discuss Hawking radiation when the mag- 
netic field is absent. In Figs. [TU1 QTJ [TJl one can see the 
emission rates for mass, angular momentum and charge 
per unit frequency per unit time, and, in the boxes, the 
results of integration over frequency u>, that is the total 
emission rates. When B = 0, the energy, angular momen- 
tum and charge emission rates of positively charged par- 
ticles are larger than those of negatively charged ones for 
all values of u), and, consequently, the total emission rates 
for positive particles are larger as well. When increasing 
the black hole charge Q, the gap between the positive and 
negative particles emission rates increases. Electrostatic 
repulsion of positive particles by the black hole (being 
proportional to the charge Q) enhances the emission of 



more positive particles. The total energy emission rate 
decreases as Q is growing, the same being true for the 
momentum emission rate. The total emission rates in- 
clude summation of both positive and negative particles, 
so that the most interesting correlation occurs for the 
charge emission rate: when Q grows, the charge emission 
rate, unlike the energy and momentum rates, increases. 
In general in geometrical units, the black hole looses its 
mass more quickly than its charge, reaching thereby the 
extremal Kerr-Newman state. 

When one turns on the magnetic field, the picture of 
Hawking radiation changes drastically. First, at rela- 
tively small values of Q, the energy and momentum emis- 
sion rates of positive particles are no longer than those of 
negative ones for all uj. At some values of oj, the energy 
and momentum emission rates due to positive particles 
are smaller than those of negative ones (see Figs. [TU1 [TTJ 
[T2")) . For large values of charge Q, the gap between emis- 
sion rates of positive and negative particles increases, so 
that the intensity of emission due to positive particles 
becomes dominant again. What is more important, the 
presence of the magnetic field considerably increases the 
energy and momentum emission rates and, at the same 
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time, considerably decreases the charge emission rate, up 
to changing the sign of the charge emission rate, which 
means re-charging of the black hole. This means that, in 
the presence of magnetic fields, the black hole evaporates 
much quicker and reaches the extremal state in a much 
shorter period of time. This is quite evident if one notices 
that the Faraday induction produces an additional (in- 
duced) charge —2aMB on the surface of the black hole. 
This charge is opposite to the black hole charge Q and 
attracts positively charged particles and repulses nega- 
tively charged ones. At sufficiently large values of the 
magnetic field, the absorbtion of positive particles will 
dominate over the negative ones, which leads to increas- 
ing instead of decreasing the black hole charge during the 
evaporation process. This process considerably decreases 
the time needed by a black hole to reach its extremal 
state. This may be a relatively small effect for astro- 
physical black holes but is not negligible for miniature 
black holes. 

Finally, let us recall that, if the induced electric field is 
as strong as /K 2 /e, the electrodynamic Schwinger mech- 
anism of pair creation will occur. Unlike the Hawking 
radiation which occurs on the edge of the black hole, the 
Schwinger process will be contributing in the particle pro- 
duction outside the black hole horizon. The Schwinger 
production will make positive particles move toward the 
black hole horizon and make negative particles move out- 
wards. Although we have not done any estimates for this 
process, qualitatively the Schwinger production should 
probably enhance the recharging of the black hole and, 
in this way, make the evaporation even quicker. 

VI. CONCLUSIONS 

We have considered the quasinormal modes, classical 
scattering (through calculations of reflection/transition 
coefficients) and Hawking radiation of Kerr-Newman 
black holes immersed in a homogeneous magnetic field. 
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